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Abstract 

Di-gluon fusion mechanism might account for the large branching ratios of 
B — * Kn' . But because we know little about the effective rj'gg vertex, there 
are large uncertainties in perturbative QCD estimations. In this paper, we 
try several kinds of rj'gg form factors and compare the numerical results with 
the experiment. We find that, though we know little about rj'gg form factor, 
if di-gluon fusion mechanism is important in B —* Krj' , the branching ratios 
of the decays B — > K*rj' would be around 10 -5 which can be tested by future 
experiments. 
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I. INTRODUCTION 



Recently, CLEO collaboration have improved their previous measurements of B — » Krf 
0: ' 

B(B+ -> = (80^° ± 8) x KT 6 

£(5° -> = (88±i| ± 9) x KT 6 . (1) 

These branching ratios are much larger than the estimations under the standard theorectical 
frame which is based on the effective Hamiltonian and general factorization approximation. 
Now it is commonly believed that these large branching ratios are due to the special prop- 
erties of r)', and several new mechanisms have been proposed to enhance the decay rates 
of B — > rj'K and B — > r)'X s . In the following we only discuss the standard model(SM) 
mechanisms because we think that the contribution of the SM should be carefully examined 
first. 

Halperin and Zhitnisky proposed an interesting mechanism: rf can be directly pro- 
duced through b — > ccs — > rj's if (7/ |c7 M 7sc|0) = —if°,P% 7^ 0. But if this mechanism 
dominates, one easily find that 

B(B° -> V 'K*°) \(i\(cc) v _ A \0)(K*°\(sb) v ^ A \B )\ 2 



B(B° - r/K°) \(rf\(cc)v-A\0)(K°\(sb)v-A\B )\ 2 

\2fc,M K *A* K *(M 2 v ,)(e K * ■ P B )\ 2 (Al 



BK* 




m{Mi- M2 M BK m)\ 2 ~\Fo BK m) < -°- 9 ' (2) 

which is in contradiction with the stringent upper limit of B° — > rj'K* given by CLEO |]J: 

B(B° -> rj'K* ) < 2.0 x 1(T 5 . (3) 

So this mechanism is unlikely to be dominant. Because the effective vertex of b — > sgg in 
SM is very small, b — > sgg srj' is also impossible to account for the large B — > r(K 
branching ratios. The authors of Ref [|J proposed that b — ► sg* and g* — ► ^77' via QCD 
anomaly can account for the large semi-inclusive branching ratios of B — >• CLEO 

have measured the 77' momentum spectrum in i? — > r]'X s and this measurement favors the 
mechanism of Ref M. But because it has an extra gluon in the final state, unless the gluon is 



soft and absorbed into the rf wave functions fll2 |, it seems difficult to realize its contribution 
to two-body exclusive decay B — > rj'K. 

Motivated by the idea of Ref [§,0, the authors of Ref [0 proposed di-gluon fusion mech- 
anism, which is depicted in Fig.l. Because both of the gluons are hard, it seems reasonable 
to give a perturbative QCD estimation. But because of our ignorance about the form factor 
of rj'gg vertex, there are large uncertainties in calculations. So in the following we would 
reanalyze di-gluon fusion mechanism in detail using several kinds of (gg\rj') form factors and 
compare our numerical results with experimental data. 



II. ANALYSIS OF DI-GLUON FUSION MECHANISM 

We first give a brief introduction on how to estimate contributions of di-gluon fusion 
mechanism in perturbative QCD. 
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g — g — rf vertex can be parameterized as || : 

(g a g b W) = 45 ab g 2 s A v ,F(kl k^e^k^e", (4) 

where A v > is a constant and can be extracted from J/\P — * 77/ or 7/ — > 77, g s is QCD 
coupling constant, F{k\, k 2 ) is (gg\r]') form factor and -F(0, 0) = 1. Neglecting the transverse 
momentum of quarks and taking the wave functions for B and K meson as 

1 I c 1 I c 

= 7/=--/= <t>B{x)U> + m B )j 5 , V K = -j=-j=(f) K (y)^ 5 , (5) 

where Ic is an identity in color space, and the integration of the distribution amplitude is 
related to the meson decay constant, 

(f) B (x)dx = — J (j) K {y)dy = — y=. (6) 



In this paper, we take the distribution amplitudes as [IC 



(^Bix) = ^=5(x -e B ), (j) K {y) = \J~y(l -y)- (7) 
Then we can write the amplitude of Fig.l as 

M= J dxdy<p B {x) e ^/Wk? 

where is the effective b — > sg vertex || 

^ = ^^- 2 V*V lb st a [Fi{kh, - - FfarKnJ-^-} 6, (9) 

and the color factor C F = -^-^Tr{t a t b }5 ab = §. Finally we obtain 
M = < =^=a 2 s A ri ,C F 32V*V tb J dxdy<f) B (^K(y)F(k 2 1 ,k 2 2 ) 



F{k 2 [{p ■ kjjq ■ h) - (p ■ k 2 )(q ■ h)} + F^m B m b [(q ■ k 2 )k\ - (q ■ kjjh ■ k 2 )] 
X Wk 2 ' [ } 



In the above integral, there are two singularities from the gluon propagators and sometimes 
there would be another singularity from the (gg\i]') form factor (For example, F(k 2 ,k%) = 

Til?' 

2(fci'fc2) )' ^ e ac ^ a sma ^ P ure imaginary number ie (e > 0) to the denominator to get a 
convergent integral. 

To evaluate the numerical result of Eq.(10), we take the form factors Fi and F 2 according 
to ref HJI3 as 

F x = — (C 4 + C 6 ), F 2 = -2C 8 . (11) 
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Where Cj are wilson coefficients at the NLL level. (Because QCD corrections maybe large, 
we do not take |||8| F\ ~ —5 and F 2 — 0.2 which are derived from the SM without QCD 
corrections.) It is always subtle to choose the scale of a s , because the average momenta 
squared of the gluons are 

(k 2 ) ~ 12 GeV 2 (kj) ~ 1 GeV 2 (12) 



we prefer to choose a s = ya^kjja^k^} = 0.28 though we also give the branching ratios 
when a s = oi s {k\) = 0.21 or a s = a^/c^) = 0.38 in Tab.l and Tab. 2. 

As to the (gg\rj') form factor F(k 2 ,k%), because of our complete ignorance there are large 
uncertainties in the numerical evaluation of Eq.(10). In the following we try several kinds 
of form factors. 

In Ref [d||, the authors have assumed that g 2 -dependence of the form factor is weak 
and as an approximation they take form factor as a constant to explain large semi-inclusive 
decay B — > X s r]'. The difference between Ref f|] and Ref || is that the running of a s {ji) with 
the scale is considered in || but not in [Q. We also use their ansatz to estimate exclusive 
decays B — > Krj' and B — > K*T]'. Using Eq.(10), we give numerical results in Tab.l. From 
the Table it seems that the constant form factor can account for the large branching ratios 
of B — > Krj', but unfortunately we would get much larger branching ratios of B — > K*rj', 
which is strongly in contradiction with the CLEO's measurements IJ (see, for instance, the 
case of €b = 0.06 in Tab.l). 

From Eq.(10), we notice that the amplitude must be integrated over a wide range of k\ 
and k 2 , therefore the effects of k 2 , k\ dependence of F(kf, k 2 ) must be taken into account. 

The authors of Ref choose the form factor &s F(k 2 , k%) = 2(ki'k 2 ) because such a form 



factor works well in J/^> 777' |H |. We have followed their calculations and it seems 
that their numerical results are overestimated than ours(in Tab. 2) by a factor about three. 
From Tab. 2, we can see that our estimations of BiB — > rj'K) are about 10~ 5 , but when 
considering that the standard theoretical frame (based on effective Hamiltonian and general 
factorization approximation) can give JRI BiB — > i]'K) ~ 3.6 x 10~ 5 , it is still possible to 
account for the experimental data only if the contributions from di-gluon fusion mechanism 
and the standard theoretical frame constructively interfere. If this is true, then because 
the contributions to BiB — > rj'K*) from the standard theoretical frame are negligible, the 
di-gluon fusion mechanism would be dominant in B — > rj'K*, i.e., BiB — > rj'K*) ~ 10 -5 
which can be tested by future measurements of CLEO or B factories. 

Kagan and Petrov |TT| proposed a model of the g — g — rj' vertex in which a pseudoscalar 
current is coupled to two gluons through quark loops. Their perturbative calculations yield 
a form factor: 

F(k 2 ,kl)(X a f m S 

f=u,d,s 

1 1-x 

ay 



C y m 2 — (1 — x — y)(xk 2 + yk^) — xym 2 , — ie 







where we take a u — aa — 1, a s — 2 and normalize F(k 2 ,k 2 ) with the normal condition 
F(0, 0) = 1. We use this form factor in Eq.(10) with e# = 0.06 and get 
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B(B -> rj'K) = 1.28 x 1(T 7 . (14) 

which is too small to account for the experiments. 

In Ref ||16|| , the authors calculate the transition form factor in 7r° coupling to 7*7* in the 
frame of a perturbative QCD based on the modified factorization formula. They find that 
numerically their results are extremely similar to that obtained by applying the interpolation 
procedure in the manner of Brodsky and Lepage in the case of two off-shell photons: 

F n o rr (klk 2 2 )^(l~^)-\ (15) 



7T 



with 

X 2 = (fcf - kl) z 

and A^ = 27r/ 7r = 0.83GeV '. We assume that the structure of the transition of rj' — > g*g* is 
similar to that of 7r° ^7*7*: 

n^^) = (i-5-r 1 , (17) 

where we approximate A^/ = A„.. By using this form factor, we calculate the branching 
ratios of B — > Kr/ which are listed in Tab. 3. The results are about fifty times smaller than 
the experimental data. 

In experimental fit, pole approximation is often used to fit the momentum square de- 
pendence of form factor [|15|]. As a try, we also assume F(fc^,fc|) = j^z^rp^r^z^rp^i ) for 



77 * 77 

our calculations, we list the results in Tab. 4. And we can see that this kind of form factor 



would make the digluon fusion mechanism completely negligible in B — > rj'K. 

To examine the validity of perturbative QCD in the above caculations, or in other words, 
whether the amplitude calculated by perturbative QCD is dominated by hard gluon con- 
tributions, we set cut on k\ and compare the results with those without cut. Taking 
Aqcd — 200MeV, we list the numerical results in Tab.V. We can see that hard gluon 
contributions are really dominant in the case of constant form factor. But if we take form 

2 

factor as F(k\,k%) = 2 ki-'k 2 ' narc ^ gluon contributions are small, this does not mean that 
non-perturbative contributions are dominant in this case. This is due to the fact that the 
singularity in the form factor is accidently close to the other singularities in k\ and k\ and 
then enhance the contributions of soft gluon. 



III. REMARKS AND CONCLUSIONS 

The authors of Ref proposed digluon fusion mechanism to explain the large branching 
ratios of B — > rj'K, but because of our ignorance about effective r]'gg vertex, there are large 
uncertainties in perturbative QCD estimations. We try several kinds of (gg\rj') form factors, 
and through our calculations, we find that constant form factor is consistent with the data 
of B — > rj'K, but inconsistent with the data of B — > rj'K* in some cases of different a s or e# 
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. if we take F(kf, fcf ) = 2fc v ' k as the authors of Ref @ have done, di-gluon fusion mechanism 
is important in B — > 77' A" but not dominant. As a consequence, we could anticipate that 
B(B — » r/A*) would be about 1CT 5 which can be tested by future experiments. We also 
try the other three kinds of form factors and they all give very small contributions to the 
branching ratios. 

We conclude that, though we know little about (gg\r]') form factor, if di-gluon fusion 
mechanism is important in B —>■ rj'K, the branching ratios of B — > r]'K* would be definitely 
around 10 -5 . 
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TABLES 





0.05 


0.06 


0.07 


Exp. § 


a s 


0.21 


0.28 


0.38 


0.21 


0.28 


0.38 


0.21 


0.28 


0.38 


BU{B- -> K~r/) 


1.99 


6.60 


21.3 


1.32 


4.40 


14.2 


0.94 


3.13 


10.1 


8.0^ ±0.8 


BK(B° -► if V) 


2.06 


6.88 


22.2 


1.38 


4.59 


14.8 


0.98 


3.26 


10.5 


8.8t[i ±0.9 




4.04 


13.5 


43.4 


2.40 


8.00 


25.8 


1.59 


5.30 


17.1 


< 8.7 


B7e(5° -» #*V) 


4.19 


14.0 


45.1 


2.49 


8.31 


26.8 


1.66 


5.52 


17.8 


< 2.0 



TABLE I. The branching ratios of B — > Krf and B — > K*r}' in unit of 10 5 by using con- 
stant f orm factor F( fc?,fc|) = 1. Where a s {k\ = 12GeV 2 ) = 0.21, a,Qs% = IGeV 2 ) = 0.38 and 
a s = ^ a s {k 2 )a s {k 2 ) = 0.28. 





0.05 


0.06 


0.07 


Exp. |] 


a s 


0.21 


0.28 


0.38 


0.21 


0.28 


0.38 


0.21 


0.28 


0.38 


BK(B- -> K-r]') 


0.60 


2.07 


6.42 


0.47 


1.63 


5.05 


0.39 


1.34 


4.15 


8.0l^±0.8 


BK(B° -> X°t/) 


0.62 


2.15 


6.67 


0.49 


1.70 


5.27 


0.40 


1.39 


4.31 


8.8±\% ± 0.9 


Bft(.B- -» K*~rf) 


0.91 


3.17 


9.83 


0.67 


2.31 


7.16 


0.54 


1.89 


5.86 


< 8.7 


BK(B° -» #*°t/) 


0.95 


3.30 


10.2 


0.69 


2.40 


7.44 


0.57 


1.97 


6.11 


< 2.0 



TABLE II. The branching ratios of B — ► Krf and B — ► K*rf in unit of 10 5 by using the 
form factor F(k\,k\) = 2{k ^' k2) • Where a s (kj = 12GeV 2 ) = 0.21, a s (k% = IGeV 2 ) = 0.38 and 
a s = y ' a s {kl)a s {k 2 ) = 0.28. 





0.04 


0.05 


0.06 


0.07 


0.08 


B(B- -► K~rf) 


2.31 


1.71 


1.34 


1.09 


0.92 



TABLE III. The branching ratios of B — ► ifr/ in unit of 10 6 by using the form factor 
F(k 2 ,k 2 ) = (l-X 2 /A 2 )- 1 . 
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0.04 


0.05 


0.06 


0.07 


0.08 


BK(B~ -► K~rf) 


2.01 


1.30 


0.91 


0.67 


0.52 


B1Z{B° -► i?V) 


2.09 


1.35 


0.95 


0.70 


0.53 



TABLE IV. The branching ratios of B -> Kr/ in unit of 10" 6 by using the form factor 

F ( k l k 2) = (l-k-(/mJ(l-k'i/m v ,)- 



Decay Modes 




Form Factor 


F(k 1 ,k 2 ) - 2(fcl ? fc2) 


F(fc 2 ,fc 2 2 ) = l 


F{k l ,k 2 ) - 2(kl . k2) 


F(fc 2 ,fc 2 ) = l 


A 2 

A QOD 


71% 


100% 


72% 


94% 




45% 


93% 


43% 


84% 


^QCD 


22% 


71% 


17% 


58% 



TABLE V. When ee = 0.06, the ratio of 'hard contribution' to di-gluon fusion mechanism 
with k\ cut of A 2 QCD , Ak 2 QCD and 9Aq CD . 
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FIG. 1. The diagram for di-gluon fusion mechanism. 
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